Multiparticle covariance matrices and the impossibility of detecting graph state 
entanglement with two-particle correlations 
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We present a criterion for multiparticle entanglement based on covariance matrices. On the one 
hand, the criterion allows to detect bound entangled states which are not detected by other criteria; 
on the other hand, some strongly entangled pure states such as the GHZ states are not detected. We 
show, however, that this is a general phenomenon: No separability criterion based on two-particle 
correlations can recognize the entanglement in the family of graph states, to which the GHZ states 
belong. 

PACS numbers: 03.67.-a, 03.65.Ud 



I. INTRODUCTION 

The presence of quantum correlations in physical states 
gives several advantages in performing certain tasks. For 
instance, entangled states are used in quantum key dis- 
tribution protocols, serve as a resource in measurement 
based quantum computation and provide insights into 
fundamental questions about quantum mechanics [l|, [^l ■ 

Consequently, many efforts have been undertaken to 
characterize entanglement. As a first step, one usually 
tries to prove inseparability of a given state. In case the 
state is entangled, one may further attempt to estimate 
the amount of the entanglement in the state. For bipar- 
tite systems, entanglement criteria based on covariance 
matrices have turned out to be powerful tools for detect- 
ing and quantifying entanglement [^-Q. 

In this paper, we present an entanglement criterion 
based on covariance matrices for multipartite entan- 
glement. Our criterion is a generalization of the co- 
variance matrix criterion (CMC) derived in Ref. Q to 
the multipartite setting. We show that the multipar- 
tite CMC is in some cases a strong criterion, detect- 
ing bound entanglement in thermal states of spin mod- 
els, which is not detected by other commonly used cri- 
teria. Surprisingly, some pure highly entangled states 
like Greenberger-Horne-Zeilinger (GHZ) states are not 
detected by the multipartite CMC. However, this turns 
out to be a quite general phenomenon coming from the lo- 
cality of the observables used in the covariance matrices. 
In particular, we show that all two-particle correlations 
of graph states (a family of highly entangled multi-qubit 
states that is important for many applications 0) are 
compatible with a fully separable state. Consequently, 
no separability criterion, which uses only the informa- 
tion gained from the two-party reduced density matri- 
ces can detect graph state entanglement. This proves 
that not only the multipartite CMC cannot detect graph 



states, but also optimal spin squeezing inequalities [8[, 
entanglement witnesses based on two-particle Hamilto- 
nians @ or structure factors or criteria based on 
magnetic susceptibility measurements [ll| will fail to do 
that. A previously noted result in this direction is that 
the entanglement of linear cluster states, which consti- 
tute a subset of graph states, cannot be detected by spin 
squeezing inequalities [l^ . 

Our paper is organized as follows. In Sec. |TT] we will 
first recall the basic notions of covariance matrices and 
of the bipartite CMC. Then, we present the multipartite 
generalization of the CMC and discuss its evaluation. For 
multi-qubit systems, one can evaluate the criterion via 
scmidefinitc programming and for tripartite systems we 
present an analytical approach. In Sec. HI we consider 
examples of states that are detected by the multipartite 
CMC. We first consider pure three qubit states, and then 
bound entangled thermal states in spin models. It turns 
out that of these states the multipartite CMC detects 
a larger subset than the optimal spin squeezing inequali- 
ties. In Sec. IV we show that the two-particle correlations 
in any connected graph state with more than two qubits 
arc indistinguishable from the correlations of some fully 
separable state. Finally, we conclude and discuss further 
research directions. 



II. GENERAL CRITERION AND ITS 
EVALUATION 



We begin this section by recalling the definition of the 
covariance matrix (CM) for the bipartite case and stat- 
ing the covariance matrix criterion (CMC) for two par- 
ties. This will serve as a fundament for the generalization 
on three parties, which can be easily generalized to any 
number of parties. 
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A. The bipartite CMC 

Let us start with the definition of CMs and some basic 
facts about the CMC. A detailed description can be found 
in Ref. Q. 

Definition 1. Let g be some quantum state acting on 
the Hubert space "H and let {Af^} be a set of observables. 
Then, the covariance matrix 7 of the state g is a real and 
symmetric matrix, given by the entries 

Ikiig) - ^{MkMi + MiMk)e - {Mk)g{Mi),. (1) 

fn what follows we will write jki instead of 'jkiig), since 
it is usually clear which quantum state is considered. 

The CMC for the bipartite case formulated in Rcfs. [1, 
Q makes use of a block form of the CM. The block struc- 
ture of the CM arises naturally if one considers only local 
observables: 



One can extend this criterion, by looking for the largest 
t, such that 7a_b(pa_b) > tKA ffl kb holds. Then, for sep- 
arable states this inequality is fulfilled for t = 1, while 
for entangled states it can happen that only smaller val- 
ues are allowed. This parameter t can then be used for 
a quantification of entanglement by giving lower bounds 
on the concurrence, for details see Ref. p. 

A central difficulty in this criterion is the characteriza- 
tion of the possible ka,b- How is it possible to guarantee 
that no KA,B exist? For this problem, however, several 
results have been obtained. Consequently, the CMC has 
been shown to be a strong criterion for bipartite entangle- 
ment: It detects many bound entangled states, contains 
several known criteria as corollaries, and, with the help 
of local filteringoperations, it is necessary and sufficient 
for two qubits [j, Q . 



B. The multipartite CMC 



Definition 2. Let g = gAB be a bipartite quantum state 
defined on the Hilbert space Ha ^ Hb of dimension d ~ 
dAdB and {A^ (8) 1, 1 ® -B/} a set of local observables. 
Then the covariance matrix of the state gAB has a block 
structure 



lAB 



A C 



(2) 



where A and B are CMs of the reduced states and the ma- 
trix Cm ~ {Ak®Bi) — {Ak){Bi) contains the correlations 
between the subsystems. For product states Q = gA® Qb 
the matrix C vanishes. 

Up to now, we did not specify the observables which 
are used in our approach. In the following, we will 
always assume that the observables Ak (and similarly 
Bk) are an orthonormal basis of the corresponding op- 
erator space. This means that they fulfill Tr{AiAj) = 
Sij] for example, if Alice has a qubit, one can take 
{Ai} = {l/^/2, ax I cry/\/2, 0-2/72}. The later results 
are, however, independent of which basis of observables 
has been chosen. 

Let us now formulate thegeneral CMC, which has al- 
ready been derived in Ref. [J| . This will make the gener- 
alization on the multipartite case more clear. 

Proposition 3 (Bipartite CMC). Let jab be a block 
CM of a bipartite separable state gAB ■ Then there exist 
pure states \ipk) in Ha and \(j)k) inH.B and probabilities 
Pk such that if we define ka = J2kPkli\'^k){'4'k\) and 
KB = Z]fcPfe7(l^fc>('/'fc|) ^he inequality 



lAB^gAB) > KA® KB 



A 

C^ 



> 



Ka 






Kb 



(3) 



holds. This means that the difference between left and 
right hand side must be positive semidefinite. If there are 
no such KA.B then the .state p must be entangled. 



Let us now start with the generalization of the CMC 
to multipartite systems. We formulate our conditions for 
tripartite systems, but the generalization to more parties 
is then straightforward. First, in analogy to Definition [2] 
we consider a CM of a tripartite state gABC by using only 
local observables. By doing so we arrive at the definition 
of the block CM for tripartite quantum states: 

Definition 4. Let g = gABC be a tripartite quantum 
state defined on the Hilbert space Ha ® Hb ® He and 
{Ak (8)l(g)l,l(g)i?;(g)l,l(8)li8) C„i} a set of local ob- 
servables. Then the covariance matrix of the .state gABC 
has a block structure 



A D E 
B F 

E^ F^ C 



(4) 



where A, B, and C are CMs of the reduced single parti- 
cle states and the matrices Dki = {Ak ® Bi) — {Ak){Bi) 
etc. contain the correlations between the subsystems. 

Similar to the bipartite case, one can easily see that for 
any fully product state g = gA Qb 9c the matrix 7 
takes a block-diagonal form. Hence, using the concavity 
property of the CMs we arrive at the generalization of 
the CMC on the tripartite case: 

Tiieorem 5 (Tripartite CM). Let g be a fully separa- 
ble tripartite state. Then there exist states l^i'''^) and 

matrices Ka of the form Ka = '}2kPkl{\'t'a^){'t'a^ \) for 
a = A, B ,C such that 



(5) 




// no such matrices Ka exist, the state is entangled. 
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Proof. The proof is similar than for the bipartite CMC 
Any fully separable state can be written as 

k 

where the {pk} form a probability distribution. Using 
now the concavity property of the CMs, stating that 
liJ^kPkS^'''^) ^ J^kPklio^'''') and the fact that for the 
product states 7 is block diagonal proves the claim. ■ 

While the multipartite CMC is easy to prove, the com- 
plicated part is again its evaluation: How can we exclude 
that the Kq, exist? In the next paragraphs we present 
two strategies to do that. Firstly, for multi-qubit sys- 
tems, one can resort to semidefinite programs (SDPs) 
and formulate the problem of searching over all possible 
/s^Q as a feasibility problem [l^ . Secondly, one can deduce 
computeable criteria from the positivity semidefiniteness 
of the matrix in Eq. ([5]). 

C. Evaluation of the multipartite CMC for qubits 
using SDP 

For qubits, the positivity test in Eq. ^ can be for- 
mulated in a simple way as a special instance of an ef- 
ficiently solvable SDP, namely feasibility problem. 
This is true for any number of qubits, but cannot be 
easily extended to higher dimensional systems. 

In order to do so we first note that in the case of 
qubits one can give an exact characterization of k^. 
Let us assume that the used observables are {Ai,i = 
1,...,4} = {l/V2,ajV2,ay/V2,<7,/V2}. Then, as one 
can directly check, the first row and column (correspond- 
ing to Ai = I/V2) of the single-qubit CM vanishes. 
So the criteria can effectively be formulated with 3x3- 
matrices k. These matrices were completely characterized 
in Ref. [H in terms of vectors in IR'^: 

Lemma 6. For any vector \(f>) g IR'^ a matrix of the form 
(I3 — |0)((/)|)/2, is a valid CM of some pure single qubit 
state. Consequently, if Alice owns a qubit the set of valid 
KA is given by all matrices of the form 

KA = ^{'i-3- Pa), (7) 

where pA is a real 3x3 matrix with trace one and positive 
eigenvalues. 

Moreover the characterization of k for one qubit pro- 
vided by the last Lemma is exhaustive in the following 
sense: 

Lemma 7. Consider the CMC condition for three qubits, 
7 ^ ka (B Kb ® Kc- Then, there exist positive matrices 
Xa,Xb and Xq such that Tr{XA/B/c) — 1 ^''^d 7 > 
Xa © Xb ® Xc if and only if there exist valid k-a/b/c 
that fulfill the CMC condition. 



Proof: The sufficiency is given by the last Lemma; 
clearly, the ka/b/c can be directly used as Xa/b/c- 
the other hand, suppose we have found such matrices 
Xa/b/c- We know, that the k^s for qubits must have 
a special form, namely = i (I3 — pi). Since the Xi 
satisfy the condition 7 > Xa © Xb © Xc and 7 has 
eigenvalues less or equal than i [H, the Xi have also 
eigenvalues smaller than ^ , and any Xi can be written in 
the form Xi ~ i (I3 — Xi), where Tr{xi) = 1 and Xi > 0. 
Therefore, Xi is exactly of the form of Eq. (I7|). ■ 

Using the last two Lemmata we can formulate the 
CMC for the three-qubit case in the following form: 

maximize: t 

subject to: 7 - tXA & Xb ® Xc >0, 
Xa,b,c > 0, 

Tr{XAM,c) = 1. (8) 

If the maximal t is larger or equal one, then the state is 
not detected by the CMC, but if the above maximization 
has a solution only for < < 1 then the state is detected 
by the CMC and is therefore entangled. 

The problem in Eq. ([5]) can be readily re-written as 
a semidefinite program (SDP). Such a problem can di- 
rectly be solved using packages like SEDUMI or YALMIP 
[l3l |. The above formulation is a simplification of a sim- 
ilar formulation for two qubits derived in Ref. Note 
that the solution of Eq. ([5]) directly allows to determine 
the parameter t, so the result might be used for further 
quantification. Finally, note that the SDP can be di- 
rectly extended to test the multipartite CMC for more 
than three qubits. Furthermore, if higher dimensional 
systems are considered, one can also test the CMC via 
a SDP (since the Hi are always nonnegative and have a 
fixed trace Q), however, this is not an exhaustive char- 
acterization of the Ki and the SDP is therefore weaker 
than the original CMC. 

The multipartite criterion evaluated by the SDP de- 
tects many interesting bound entangled states, which ex- 
amples will be discussed in the next section. For the time 
being we discuss a method of analytical evaluation of the 
tripartite CMC. 

D. Analytic test of non-separability of tripartite 
states 

In contrast to the previous approach, the analytical 
method works only for the tripartite case, but it is not 
restricted to qubits. 

What are the possible corollaries of the non-negativity 
of the block matrix in Eq. ((S])? First of all, we know that 
all submatrices of this matrix have to be non-negative. 
As we learn from the bipartite case the strategy of 
considering submatrices leads to a quite strong and easily 
computable separability criterion. In order to investigate 
whether a tripartite state is fully separable or not, it is 
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then natural to consider 3x3 submatrices of the CM in 
the block form. To do this, we need the following fact 
which is true for any real 3x3 non-negative matrix: 



Lemma 8. Let rj be a 3 x 3 real positive matrix 




(9) 



then 



a \d\ \e\ 
V = \ \d\ b I/I I >0 
|e| I/I c 



(10) 



is also a positive matrix. 

This Lemma is an example of a norm compression in- 
equality and the proof is given in the Appendix. 

Using this property of 3 x 3 non-negative matrices we 
can prove a proposition, which gives us an easily com- 
putable separability test for tripartite states: 

Proposition 9. // the state g on a d x d x d-system is 
separable, then 

(Tr{A)-d+l EJ4,| EJe»l \ 

Tr{B)-d+l Y.^\M > 0. 
V EJe.,| EJ/«I Tr{C)-d+ll 

(11) 

Here, dij,eij, ftj, are the elements of the blocks D,E,F 
in Eq. respectively. 

Proof. Since j — ka(x) kb (X" kc is nonnegative, all 3 x 3 
submatrices are positive. So we consider the matrices 

dii hu - fit I > 0. (12) 

^ii fa ^ii {^c)i 

Applying now Lemma [5] and summing over all i gives the 
claim, if we use in addition the fact that Tr{Ka) = d — 1 
for a — A,B, C, shown in Ref. [sj. ■ 



III. EXAMPLES 

In this section we consider some examples for the pre- 
viously derived criteria. We will first consider random 
pure three-qubit states, and then thermal states in spin 
models, which may be bound entangled. 



be brought by local unitary transformations into a gen- 
eralized Schmidt form, 

IV') = AolOOO) + Aie^-^llOO) + A2IIOI) + A3IIIO) + A4IIII): 
with Ai > 0, < (/) < TT, ^ A,2 = 1. (13) 

i 

Varying the real parameters A^ and (j) one obtains differ- 
ent families of entangled states. 

Proposition [H] detects almost all of these entangled 
three-qubit states: Choosing the Pauli matrices as ob- 
servables and the parameters A^ is uniformly distributed 
between and 1, whereas 9 uniformly distributed be- 
tween and TT, one is able to detect 100% of the ran- 
domly generated states. However, Proposition [9] fails to 
detect entanglement is the case of GHZ states, |GHZ) = 
(|000) + |lll))/\/2. Nontheless, even a small perturbation 
of the GHZ state avoids this problem, e.g., the state 

l^GHz) = ^ (1000) + 10-5|110) + 1111)) (14) 

is detected, here Af denotes the normalization. 

In order to test whether this is an issue of the Propo- 
sition [5] or whether is a property of the general CMC 
as stated in Theorem [5] we used the first strategy and 
evaluated the general criterion with an SDP. The SDP 
detects also 100 % of the randomly generated states (in 
both cases 10"* instances were tested), but it also fails to 
detect GHZ states. 

This phenomenon can be explained as follows: In 
the CMC as in Theorem [5] uses only one- and two- 
point correlations (i.e. expectation values like (cr^) or 
(ct^ct^)) but never higher order correlations (like three- 
point correlations (o'^o'^o-^)), as only the former appear 
in the block structure of the CM. However, consider- 
ing only one and two-point correlations, it is impossi- 
ble to distinguish | GHZ) (GHZ | from fully separable state 
( 1 000) (000 1 1 1 1 1) ( 1 1 1 1 ) /2 , as for both states the reduced 
one- and two-particle density matrices are the same. 

As we will see in Sec. II VI this is not an exceptional case. 
Any state in the large class of graph states, which also 
contains the GHZ states, cannot be distinguished from a 
fully separable state by considering one- and two-particle 
correlations. This does not only affect the multipartite 
CMC, but it also means that these states can never be 
detected by other approaches like spin squeezing inequal- 
ities [3 , or witnesses based on two-body Hamiltonians 
or structure factors [lOj . 



B. Entanglement in thermal states 



A. Random three-qubit states 

In order to test the performance of the derived en- 
tanglement criteria we first investigate three qubit pure 
states. According to Ref. [l^ any three-qubit state can 



Now we discuss thermal states of spin models and 
investigate their entanglement properties using the ap- 
proach of Eg. ([8]) and Proposition [9] 

In Ref. [15[ it has been shown that thermal states of 
spin models with a small number of spins can have inter- 
esting separability properties. This has been investigated 



5 



frequently with spin squeezing inequalities [l6| . They 
can be formulated as follows. First, one takes the collec- 
tive spin operators — \ 'Y^=\ '^nl where the vectors 
{fik}k are orthonormal. Then one defines the matrices 



Cfci — -^{JkJl + JlJk), 
^kl = Ckl — {Jk){Jl), 

X = {N-l)T + C, (15) 
Given these matrices, a separable state fulfills [l6j : 

Tr(r)>|, (16) 
N 

Amin(x)>Tr(C)--, (17) 

An.ax(x) <{N- l)Tr (r) - (18) 

where N is the number of spin-i particles. If one of 
these inequalities is violated, the state must be entangled. 
Moreover, if only the averaged two-particle correlations 
are known, then in many situations these inequalities are 
the strongest entanglement criteria possible, in the sense 
that for any state that fulfills these criteria, there is in- 
deed a separable state with the same C, T and x ■ 

The thermal states, which were considered in Refs. d, 
[Tsl . [l6j , are the thermal states of the Hamiltonian 



H ~ Si ■ S2 -\- S2 ■ S3 + Si ■ S3 



where Sk — io'x^\o'y^\cri^'')- This Hamiltonian describes 
three spin-^ particles interacting via the Heisenberg in- 
teraction and subjected to an external homogeneous 
magnetic field of the strength h. The thermal states 



(19) 



grih) 



-H/T 



Tr(c 



-H/T\ 



(20) 



represent a two-parametric family of density matrices. 
Using the spin-squeezing inequalities one can find regions 
in the T-h diagram where all thermal states have a pos- 
itive partial transposition (PPT) with respect to any bi- 
partition (and therefore do not violate the PPT criterion 
[l7|), but the spin squeezing inequalities show that the 
state is not fully separable. For h = and certain tem- 
peratures the thermal state is even separable with respect 
to any bipartition, but not fully separable [1, [l^ . 

As a first example, we can test our general criterion 
Eq. ([5]) via semidefinite programming for these thermal 
states. It turns out, that the multipartite CMC detects 
exactly the same amount of the states as the spin squeez- 
ing inequalities. This can be understood as follows: On 
the one hand when the spin-squeezing inequalities are 
applied to spin models, then usually Eq. (fTS)) detects the 
entanglement. This inequality, however, can be seen as 



derived from local uncertainty relations for special ob- 
servables [1, [ll[ . On the other hand, the bipartite CMC 
is known to be equivalent to local uncertainty relations 
with arbitrary observables Q and this connection clearly 
holds for the multipartite case. So it is not surprising 
that for this spin model, the CMC is not weaker than 
the spin squeezing inequality. 

However, slightly changing the Hamiltonian gives 
rise to yet another family of the states with two param- 
eters, for which the comparison of three criteria (PPT, 
Spin Squeezing and CMC) shows that the CMC criterion 
can detect some states that are not detected by the PPT 
criterion or the spin squeezing criteria. This change of 
the Hamiltonian consists in changing of the direction of 
one of the local magnetic fields, therefore destroying the 
spatial symmetry of the Hamiltonian 



H' = Si-S2 + S2-S3 + Si-S3 + h (al +al + al) 



(21) 



Detection of these states by the CMC via Proposition[9] 
or the SDP in Eq. ^ is presented in Fig.[TJ As expected, 
the Proposition [9] is weaker than the SDP. The detection 
of the PPT criterion and the spin-squeezing inequalities 
is presented in Fig. ^ Part (a) corresponds to the PPT 
criterion. Here we characterize states by taking the max- 
imum of three values of the negativity corresponding to 
three different bipartitions [l^- Part (b) corresponds to 
the detection of these states by the spin squeezing in- 
equalities. Comparing Fig. 2(a) and Fig. 2(b) one sees 
that spin squeezing inequalities are able to detect bound 
entanglement. More importantly, however, by compar- 
ing Fig. 1(b) with Fig. 2(b) one sees that the CMC 
detects more states than the spin squeezing inequalities, 
which one sees looking at the region of the plots where 
kT ^ 7, h ^ 12. However, the CMC detects better 
not only bound entangled state, an example of a bound 
entangled state, which is not detected by spin squeezing 
inequalities but violates the CMC is a thermal state with 
parameters kT = 5.533, h = 4.3. Note that this is not 
a contradiction to the optimality of the spin squeezing 
criterion for some cases, since the CMC does not only 
use the averaged two-point correlations. 



IV. ENTANGLEMENT IN GRAPH STATES 
AND TWO-PARTICLE CORRELATIONS 



In Section IIII Al we have seen that the multipartite 
CMC is unable to detect the entanglement in a three- 
qubit GHZ state. The reason was that the CMC in- 
volves only two-particle correlations and that the GHZ 
state (|000) + |lll))/\/2 and the fully separable state 
(|000)(000| -I- |lll)(lll|)/2 have the same reduced two- 
particle density matrices, hence the two-particle correla- 
tions are the same. 

In this section we will show that this observation is 
not a coincidence. We will consider the family of graph 
states (defined precisely below) which comprises not only 
the GHZ states, but also other states like cluster states. 



(a) 




h 

(b) 

FIG. 1: (Color online) Detection of the thermal states of the 
Hamiltonian H2ip for different values of temperature and ex- 
ternal magnetic field by the multipartite CMC. The dark red 
region corresponds to the states, which are not detected by 
the criterion, (a) The criterion of Proposition |9] with the 
Pauli matrices as observables. Different colors correspond to 
the most negative eigenvalue of the matrix in Eq. (|lip . (b) 
Evaluation with the semidefinite program in Eq. ((8|. Here, 
different colors correspond to different values of t. 

and which is of eminent importance for many applica- 
tions and experiments 0] ■ We will show that for any con- 
nected graph state of three or more qubits there is a fully 
separable state, which has the same one- and two-qubit 
reduced density matrices. Hence, graph state entangle- 
ment can never be detected by measuring two-particle 
correlations only, and this is not only a restriction to the 
multipartite CMC, but also to a variety of other entan- 
glement criteria. 

More precisely, we will prove: 

Theorem 10. Let \G) be a connected graph state with 
more than two qubits. Then, there exists a fully separable 
state qg such that for any pair of qubits {i,j} the cor- 
responding two-qubit reduced density matrices of \G){G\ 
and QG coincide, 

e^A\G){G\) = q^Aqg). (22) 

Moreover, qg has also the same one-particle reduced den- 
sity matrices, Q^{\G){G\) = Qi{gG)- 



(a) 




Q £ 4 B 10 IE 

h 



FIG. 2: (Color online) Detection of the thermal states of 
the Hamiltonian (|2ip for different values of temperature and 
external magnetic field by the PPT criterion and the spin- 
squeezing inequalities. The dark red region corresponds to 
the states, which are not detected by these criteria, (a) The 
violation of the PPT criterion characterized by the negativity 
[l^ . (b) The detection of the states by spin squeezing in- 
equalities (|16|l - (|18|l . Different colors correspond to the max- 
imal violation of the inequalities p6|) - (|18[ ). Comparing this 
with Fig. 1(b) for values kT ~ 7, /i ~ 12 we immediately 
see that the CMC for three qubits detects more states than 
the spin squeezing inequalities. These states violate, however, 
the PPT criterion as well. One finds an example of a bound 
entangled state, which is not detected by spin squeezing in- 
equalities, but is detected by the CMC for the following values 
of parameters kT — 5.533, h = 4.3. 

Here and in the following, we use the definitions 
Qi{A) = Tri\,{A) and g,j{A) = rr/\,j(A), where A is 
an operator acting on the Hilbert space of all N qubits, 
and I is the set all particle indices. In a slight abuse 
of the usual notation, we use the same definition below 
even in situations where A, and hence gi and gij, are no 
quantum states. 

The point of this theorem is that many entanglement 
criteria use only two-particle correlations for the entan- 
glement detection. Consequently, these criteria must fail 
to recognize the entanglement in graph states: 

Corollary 11. Consider a graph state of three or more 
qubits. Then, the entanglement cannot be detected using 
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the following criteria: 

(a) the multipartite CMC from Eq. (0), 

(b) the optimal spin squeezing inequalities 0/, 

(c) entanglement witnesses based on two-body Hamilto- 
nians Q/, and 

(d) entanglement witnesses based on structure factors 
flOl ] or inequalities using the magnetic susceptibility fU l. 

Before proving the Theorem, let us give a short defini- 
tion of graph states, for a broader discussion see Ref. 0- 
Consider a graph with N vertices and some edges con- 
necting them. We are only interested in connected graphs 
which do not separate into two unconnected parts. For 
any vertex i we can write down an operator 

g, = X, (g) Zj, (23) 

where here and in the following, Xi,Yi, and Zi denote the 
Pauli matrices acting on the i-th qubit and Af{i) is the 
neighborhood of i, that is, all vertices that are connected 
with i. The graph state \G) is now defined as the unique 
iV-qubit eigenstate of aU the gi, 

\G)=g^\G). (24) 

The graph state is not only an eigenstate of the gi, but 
also of their products. All products of the gi form a com- 
mutative group with 2^ elements, the so-called stabilizer 
S{G). Using the elements of the stabilizer, we can also 
write the graph state as 

= ^ E (25) 

This relation will be useful in our proof, which can be 
split into three parts: Firstly, we will distinguish the 
different cases where Qij is not the identity. Secondly, 
we will investigate in detail what kind of reduced states 
gij{\G){G\) can occur. Finally, we will see that there ex- 
ists a measurement of the type Ai = O^'^ with 
O*-'-* G {Xi,Yi,Zi} which allows to compute all two- 
particle correlations of the state |G) and consequently 
all gi{\G){G\) and gij{\G){G\). For such a measurement, 
however, there is always a separable state giving the same 
results as the graph state. 

Part I — Let us start by characterizing all one- and 
two-particle reduced density matrices. First, note that 
for a connected graph one always has gi{\G){G\) = 1/2. 
This can be seen from Eq. ([25|l . For most of the Sk €! 
S{G) (which are products of the gi) we have gi{Sk) = 
since they act non-trivially on qubits besides i. Since 
Sq = 1 ^ S{G) we have for this element that gi{So) = 
2-/V-1J Furthermore, we can only have gi{Sk) 7^ and 
Qi{Sk) 7^ 2^^^1 if there is a correlation operator with 
the property gi = Xi. This, however, cannot occur in a 
connected graph. 

Let us now turn to the reduced two-qubit density ma- 
trices. As above, we find for most of the stabilizing opera- 
tors that gij{Sk) = 0. Consequently for most of the qubit 




FIG. 3: (Color online) The three different cases where the 
reduced two-qubit density matrix differs from the identity. 
In case (a) the qubit i is connected to the qubit j only, and 
the qubit j is connected to the rest of the graph R. In case 
(b) the qubits i and j are disconnected but have the same 
neighborhood Mii) = In case (c) the qubits i and j are 

connected and have the same neighborhood apart from 
i.e., M — A/'(j) \ {j} = A/'(j) \ {i}. Small dots denote possible 
multiple connections. See the text for further details. 

pairs {ij} we have that gij{\G)((G\) is completely mixed, 
i.e., ftj(lG)(G|) = 1/4. However, as already discussed 
in Ref. [I^, there are three exceptions, see also Fig. [3l 
These exceptions occur if some Sk acts non-trivially only 
on two qubits i and j. For these pairs {ij} there are three 
possibilities: 

(a) Qubit i has only one neighbor j [A/'(i) = j], which 
is coupled to the rest of the graph. Note that there is no 
symmetry in interchanging the qubits i and j. Then we 
have 

ftj-(|G)(G|) = i(l + X,Z,). (26) 

(b) The qubits i and j are not connected but have the 
same neighborhood \M(i') = A/'(j)], which is connected 
to the rest of the graph. This case is symmetric in i and 
j. The reduced state is 

Q^A\G){G\)^\{\^X,X,). (27) 

(c) The qubits i and j are connected and have the 
same neighborhood otherwise which is connected to the 
rest of the graph. Defining A/'(i) ~ A/'(i)\{j} and ^^(j) = 
A/'(j)\{i} this means that M{i) = -^(j)- This case is also 
permutation symmetric. We arrive at 

ftj(|G)(G|)==i(l + y,y,). (28) 

Part II — Let us now discuss the relations between the 
three cases above. Especially, we would like to investigate 
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in which cases one particle can contribute to two or more 
cases. 

(i) For a given pair of qubits {io, Jo} maximally one of 
the three cases can apply. 

Proof. Clearly, the case (b) excludes the case (c) 
and vice versa, since in (c) the qubits are connected 
and in (b) not. The same holds for (a) and (b). 
If (a) and (c) hold simultaneously, then = 
Af{j) = which means that the graph consist only 
of two qubits or is not connected. Both cases were 
excluded. 

(ii) If the case (c) holds for some pair {ic,jc}, then 
neither ic nor jc can be involved in a case (a) . 
Proof. Assume ic = ia- Then, since ia has only one 
neighbor we have M{ic) = = M{ic) which is not 
allowed [see also (i)]. If ic = ja then ia G M{ic)- 
We cannot have ia = jc [see (i)] but we can also not 
have ia S Af{jc), since ia has only one connection. 
This is then a contradiction to J\f{ic) = ■N'ijc)- 

(iii) If the case (c) holds for some pair {ic,jc\^ then nei- 
ther ic nor jc can be involved in a case (b). 
Proof. Assume ic = ib, then jc =/= jb since prop- 
erty (i) holds. We have that jb 4- ■^i'ib) = A/'(ic) 
which implies that jb ^ -^(jc)- On the other hand, 

we have that jc € A/'(ic) = ^fiib) = J^ijb)- This 
means that both ib and jb are connected to jc, hence 
jb G -^{jc), which is a contradiction. 

(iv) If the case (b) holds for some pair {ib,jb} then there 
exists no pair {ia,ja}, for which the case (a) holds, 
such that ib = ja- But there might exist a pair 
{*a;io}) of case (a) such that 4 ~ ia- 

Proof. Assume that if, ~ ja. Then ia G J\f{ib), but 
ia € J^{jb) cannot be true since ia has only one 
connection to the rest of the graph. So M{ib) 7^ 
M{jb) which is in contradiction to the case (b). The 
second part of the property is proved by providing 
an example. For this, one can just take a linear 
graph with three vertices, numerated from left to 
right. Then, the first and the second qubits are ia 
and ja respectively, whereas the first and the third 
qubits are ib and jb- 

(v) For a pair {ia, ja} of case (a), the qubit ia may also 
appear in a case (b) . The qubit ja cannot appear in 
a case (b) or (c], but it may also appear in another 
case (a) pair {ia,ja\ but only as ja- 

Proof. The first part just follows from the previous 
points. For the last part, one may consider the 
example in (v), there, the second (middle) qubit 
takes part in two different cases (a). It is clear, 
that ja cannot be ia for some other case (a), unless 
we have a two-qubit graph. 

(vi) For a pair {ib,jb}, both if, and jb can take part 
in another pair of case (b). In analogy, for a pair 



{ic, jc}, both ic and jc can take part in another pair 
of case (c). 

Proof. We just give examples where this may hap- 
pen. For case (b), consider a graph where one cen- 
tral qubit has a single connection to to > 2 other 
qubits. For case (c) consider to > 3 qubits which 
are all interconnected. 

These are all interesting cases we have to consider. In 
summary, this proves the following: A single qubit k can 
only take part in two different cases {i,j} and {TO,n}, 
if the nontrivial Pauli matrix contributions from Qij and 
Qmn in Eqs. ((26l [27l [28|) have the same Pauli matrix on 
the qubit k. For instance, in property (iv), if fc = i^ = 
ib, then in Eqs. (^51 Wf\ both Pauli matrices on qubit j 
have are the same, namely X^. The same is true for the 
property (v), here, when k = ja ~ ja, the Pauli matrix is 
Zk. If (vi) occurs, then the observable will be Xk in case 
(b) and Yfc in case (c). 

This implies that for an arbitrary graph with more 
than three qubits, we can choose a fixed observable 
O''"'' G {ATfc, Yfc, Zk} for any qubit k, such that any non- 
vanishing two-body correlation can be represented by a 
tensor product of the O'-'^^ Note that since just a single 
observable O^'^^ occurs on each qubit, it can be changed 
to Zk with a local unitary operation, which does not in- 
fluence the entanglement properties of the graph state 

Part III — Finally, we can prove the claim of Theorem 
10. Consider the observable 

N 

Mg^(^0^''\ (29) 

k=l 

which is constructed for a special graph state \G). 

This observable has 2^ eigenvectors (which are prod- 
uct vectors) and measuring this observable results in 2^ 
probabilities, which can be labeled by the possible results 
±1 for the N qubits. In order to do this, we introduce 
a vector a where the entry ak is equal to if the re- 
sult on qubit fc is +1 and equal to 1 is the result on 
qubit k is —1. Wc denote the eigenvectors by |d;). Note 
that they are product vectors |cf) = <Sk\ctk) defined by 

OW|afc) = (-l)"H«fe)- 

With this notation, we can write the expectation value 

of the observable Ma as {Mq) = Yl,s(~^)^''°"'Ps^ 
where ps = \{d\G)\'^. Further, from the measure- 
ment probabilities ps, it is possible to compute all 
non-vanishing two-point correlations of the graph state 
\G), since they are, as mentioned above, of the type 
(0(0(90"))^ where both O^*) and O^^) appear in the ob- 
servable Mg- 

Howcver, the fully separable state gfs = '^sP3\(^){oi\ 
with exactly the same probabilities ps will reproduce 
all the non- vanishing two-body correlations of the graph 
state. In fact, it also reproduces the vanishing ones: 
Let us assume for simplicity that Aic = Z®^ . Then, 
the state can be chosen such that it will reproduce 
all correlations (also vanishing two-body correlations or 
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the non- vanishing higher order correlations) of the graph 
state in the basis defined by A^g = Z"^^ . Other two- 
body correlations, such as (XiYj), will vanish for the 
graph state. However, we can expand gfg in terms of 
tensor products of Pauli matrices as 

gfs= f3k^...kr.z'''®z'''®...®z''^ (30) 

ki=0,l 

since \a){a\ = (1/2^) {g)f^i [1 + {-l)"''Zk] can be ex- 
pressed in this form. 

This implies, that for gts the two-body correlations 
(and also higher correlations) in a different basis than 
A4g = Z®^ vanish. All in all, the fully separable state 
Qfs has the same one- and two-qubit correlations than the 
graph state, so the reduced density matrices must be the 
same. I 



V. CONCLUSION 

In conclusion, we have derived a generalization of the 
covariance matrix criterion to the multipartite case. We 
have seen that the resulting criterion is strong, but some 
highly entangled states are not detected. This is, how- 
ever, an example of a more general phenomenon, which 
is that no criterion based on two-particle correlations can 
detect graph state entanglement. 

There are several ways to extend our results in the 
future. First, one can aim to develop a criterion based 
on covariance matrices for the detection of genuine mul- 
tipartite entanglement. The present criterion can only 
exclude full separability. Second, one can try to connect 
the entanglement parameter t in the formulation of the 
multipartite CMC in Eq. ([8]) with multiparticle entangle- 
ment measures, similarly as it has been done in for the bi- 
partite CMC [g] . Finally, it would be interesting to gain 
further insight into the question which types of entan- 
glement can be verified by two-particle correlations and 
which not. This question is also experimentally relevant, 
since in some experiments only two-particle correlations 
can be measured. 



Acknowledgements 

We thank J. Eisert, B. Jungnitsch, M. Kleinmann, 
C. Knapp, S. Niekamp, N. Liitkenhaus and G. Toth 
for discussions and acknowledge support from the 
FWF (START prize and SFB FOQUS) and the EU 



(NAMEQUAM and QICS). This work was finished after 
O.C. moved to the IQC in Waterloo, and he is grateful 
for support from the Industry Canada and NSERC 
Strategic Project Grant (SPG) FREQUENCY. 



APPENDIX 

Proof of Lemma [SI' We shall split the proof into four 
parts: 

(i) If all off-diagonal elements of 77 are positive, the 
claim follows immediately. 

(ii) If exactly two off-diagonal elements of 77, say, d 
and / are negative then, we have that t) — DrjD^ with a 
diagonal matrix D = diag(l, — 1, 1), so the positivity of 
77 is equal to the positivity of fj. 

(iii) If only one element is negative, say e < 0, we have 
to prove that positivity of 77(e) guarantees the positivity 
of 77 = 7/(|e|). To this end we show an equivalence of the 
positive semidefiniteness of any 3x3 matrix A and pos- 
itive semidefiniteness of any of its 2 x 2 submatrix and 
its determinant. It is clear that the implied conditions 
on the submatrices and the determinant are necessary. 
Sufficiency is proved considering the characteristic poly- 
nomial of a 3 X 3 matrix A, which is given by 

3 

Xa(A) = -A-'^ + Tr(A)A2 - ^ det(A,,)A + det{A). 

i=l 

(A-1) 

Here, An is the 2x2 submatrix of A, which arises if the 
i-th column and row are deleted. 

Since A is symmetric, it has real eigenvalues and all 
roots of this polynomial are real. If det^Au) > and 
det{A) > 0, then the zeroes clearly have to be nonnega- 
tive, because then for A < one has xa(A) > 0. 

So let us consider 

det(77) = abc + 2e\f\\d\ - he^ - aj/p - c|dp 

= abc - 2|e||/||d| - be^ - a|/p - c\d\^ > 0. 

(A-2) 

From this it follows that dct(77) = det(77) + 4|e||/||d| > 
det(77) > 0, moreover, the determinants of 2 x 2- 
submatrices of 77 arc the same as the ones from 77 and 
hence positive. Therefore, 77 must be positive semidefi- 
nite. 

(iv) Finally, if all off-diagonal elements of rj are nega- 
tive we can first flip the sign of two of them as in (ii), 
and then we arrive at the third case. ■ 
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